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ABSTRACT 


In recent time a single comprehensive model for driver behavior has been developed using potential 
field theory of the robot motion planning. These models, at present, do not have any mechanism for 
calibrating their parameters. In this thesis an attempt is made to develop a mechanism which can be 
used to calibrate the parameters of these models. Results using the proposed mechanism show that 
the proposed mechanism can form the ground work for the a larger and more extensive calibration 
mechanism for the potential field theory based model of driver behavior. 



Chapter 1 


INTRODUCTION 


Microscopic models of traffic flow form the foundation for the theory of traffic flow. Microscopic 
models attempt to predict the drivers behavior in small interval of time. Driver behavior is captured 
in these models through to variables: (i) the steering angle chosen by driver, and (ii) the acceler- 
ation/deceleration rate use by the driver. Although, various microscopic models exist (Pipes [4], 
Forbes [5, 6], Chandler et ah [7], Herman et al. [8, 9], Gazis et al. [10, 11], Rockwell et al. [12], May 
and Killer [13], Kikuchi and Chakroborty [14], Chakroborty and Kikuchi [15], Taragin [16], Michaels 
and Gozan [17]) only two of them, namely, Vasishta [2] and Agrawal [1] attempt at developing a 
comprehensive model of driver behavior. That is, unlike the others, these models attempt to capture 
driver behavior in various driving situations (like, car-following, overtaking, free flow, etc) through 
a single model. These models are based on the potential field theory of robe/ motion planning (see 
Latombe [3]). However Vasishta [2] and Agrawal [1] do not suggest how various parameters of their 
model can be calibrated given data on observed driver behavior. This thesis attempts to overcome 
this short coming by developing a calibration algorithm. 

The thesis is divided into six chapters of which this is the first. The second chapter describes 
the problem tackled here and lays down the scope of present thesis. The third chapter presents 
results from a detail parametric study of the potential field theory based models. The fourth chapter 
describes the proposed calibration mechanism. The fifth chapter presents the results obtained using 
proposed calibration mechanism for various different scenarios. The sixth and last chapter concludes 
this thesis by summarizing the achievement of this thesis and by discussing the short comings of the 
present work. 
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Chapter 2 


PROBLEM STATEMENT AND 
SCOPE OF STUDY 


In this thesis the problem of calibrating the parameters of potential field functions of potential field 
based theory model of traffic flow is studied. However, before presenting a detail description of 
problem and scope of study, a brief description of potential field theory based model of traffic flow 
is given. 


2.1 Potential field theory based model 

Vasishta [2] and Agrawal [1] have over last three years developed the potential field theory based 
model of traffic flow. This is a microscopic model which basically describe drivers’ action in terms 
of steering angle choice and acceleration/deceleration choice under different driving scenarios. In 
this models a driving scenarios is defined as a conglomeration of traffic and roadway features in the 
immediate vicinity of a driver. Roadway feature, like road edges, curvature, potholes, etc, and traffic 
features, like parked vehicles, moving vehicles in the same and opposite direction, are all assumed 
to be obstacles which a driver avoid while driving. Each of these obstacle are assumed to emanate 
a repulsive force field around it. This force field is referred to as a potential field and is described 
mathematically through a potential field function. 

The predicted actions of a driver, from this model, depends largely on the potential field 
functions of the obstacles. Despite this depends no attempt are made either by Vasistha [2] and 
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Agrawal [1] to develop the mechanism for calibrating the parameters of these functions. This lack of 
an available mechanism formed the motivation for this study. 


2.2 Problem statement and scope of study 


The problem here, as stated earlier is to develop a mechanism which can be used to calibrate the 
parameter of potential field functions. The mechanism should take as input observed path of a 
driver in a real world situation and the predicted path from the model for the same situation with an 
initial assumed parameter values. The mechanism should produce as its output a set of calibrated 
parameter values which when used by the model produces the predicted path similar to the observed 
path. 

In order to develop such a mechanism two tasks need to be accomplished. First, a study on 
the effect of different parameters on the predicted path need to be conducted. Second, based on the 
information available from this study, rules of calibration which tell us how to change the parameter 
values need to be developed. These tasks formed the subject' matter of next two chapters. 

Given the complex nature of potential field functions and the complex manner in which they 
affect the predicted path the scope of this study is kept limited. Basically in this study the scope 
is limited to calibrating the parameters of road edges and static obstacles likes potholes, parked 
vehicles, etc. 



Chapter 3 


PARAMETRIC STUDY OF 
POTENTIAL FUNCTIONS 


As stated earlier, the tasks to be undertaken in this thesis are (i) to study the effect of potential 
field function parameter values on the driver’s position profile and (ii) to develop an algorithms for 
modifying these parameters so that the predicted path emulates the observed path better. The first 
of these two tasks is the subject of this chapter. The chapter is divided into three sections. The first 
presents the potential field function for the obstacles (or driving features) which are considered in 
this thesis. The second and third sections present the effect of the parameters of the potential field 
functions. 


3.1 Potential field functions for the different driving features 


The driving features (i.e. roadway or traffic features) considered in this thesis are: (i) road edges 
and (ii) static obstacle (like potholes, parked vehicle etc). In this section the potential field function 
for these features as used by Vasishta [2] are presented. This is done so as to make the description 
of the parametric study complete. 
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3.1.1 Potential field function for road edges 

The potential field function for a road edge is given by 

U e = ae~ hz (3.1) 

where, z is the distance from the road edge, a is a parameter which denotes the potential at the edge 
and b is a parameter which determines the rate of decrease of the potential with z. Specifically the 
function for the left edge, Ule, and the right edge, Ure, can be written as 

Ule = (3.2) 

and 

Ure = a 2 e- b ^ w ~ x) (3.3) 

where, x is the distance from the left edge of the road; a\ and a 2 are the potential values at the left 
edge and right edge respectively; bi and 62 are the slope parameters for the left edge and right edge 
potential field functions respectively. The combined potential due to both the edges of the road, Ube, 
can be written as 

U be = ai e~ hx + a 2 e- b ^ w ~^ (3.4) 

Note that it is meaningful to only work at Ube rather then then Ure and Ule since a road always 
. has two edges. Further as suggested by Vasishta [ 2 ] 


• &i > 0, & 2 > 0 

• <22 > ai ( because, in India drivers follow a “keep left” policy) 

3.1.2 Potential field function for static obstacles 

The potential field function for a static obstacle is given by 

U. = ai e-^^ 2+ ^ 2 -^ 


(3-5) 


where 


P = 



jf (£^il) 2 + (£U!0) 2 > i 
V mi ' 1 V. ma ' — 

otherwise 



6 


and 

( 64 if y < h 

f(y) = l ay 3 4* py 2 + 7y + 6 if t x < y < t 2 
[ 65 if y > t 2 

where (x,y) is the coordinate of any point on the road, (xo, yo) are the coordinates of the center of 
the obstacle; m;, m a are respectively the semi-minor and the semi-major axes of the ellipse which 
circumscribes the obstacle; p is a boolean variable representing whether point (x,y) is inside the 
ellipse (p = 0) or not (p = 1). Further, a 4 is the maximum value U Q can take and f(y) is a function 
which determines the slope for U 0 . t\ and t 2 are lower and upper limits of y with which f(y) varies. 

Before leaving this section clarifications regarding (i) the determination of the coefficients of f(y) 
and (ii) m t need to be given. 

Clarification (i) 

Given the value of 64, 65 and t 2 the coefficients a,/?, 7, <5 are calculated using the following 
restrictions: (i) f(y) should pass through (£1,64) and (^2,^5), and (ii) f(y) should pass through a 
point tp where its slope is also zero. Therefore, it needs to noted that tp is also a parameter which 
affects the predicted path. 

Clarification (ii) 

The parameter is not considered as a calibration parameter, since Vasishta [2] suggests a procedure 
for selection of m l from the size and shape of a obstacle. Hence given the obstacle, rrii is fixed. 

3.2 Effect of parameters of Ub e on the predicted driver path 

The following are identified, as indicated by Vasishta [2], as parameters of Ube • ai, a 2 , 61 or b 2 and 
The effect of these parameters on the predicted path are studied here. 

Effect of ai 

Figure 3.1 shows the effect of a\ on predicted driver path on a 6 m road keeping all other parameters 
as constant ( a 2 = 4.0, bi — 0.32641, b 2 — 0.04663 and = 7.0). It can be seen that, when the 
value of a\ increases the predicted path shifts towards right edge of the road. The same behavior 
can be observed from other road widths also. 
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Effect of a 2 

Figure 3.2 shows the effect of a 2 on predicted driver path on a 6 m road keeping all other parameters 
as constant (ai = 1.0, b\ = 0.32641, b 2 = 0.04663 and = 7.0). It can be seen that, when the value 
of a 2 increases path shifts towards left edge of the road. The same behavior can be observed from 
other road widths also. 

Effect of b\ or b 2 

Figure 3.3 shows the effect of b 2 on driver path on a 6 m road keeping ratio of and all other 
parameters as constant (ai = 1.0, a 2 = 4.0 and fj- = 7.0). It can be seen that, when the value of b 2 
increases the predicted path shifts towards left edge of the road. It should be noted that keeping 
constant and changing b 2 also implies changing bi in the same direction (i.e. increasing b 2 implies 
increasing 6i). Similar behavior can be observed for different road widths and r 1 ratio. 

t>2 

Effect of 

The effect of this can not be studied separately from effect of &i or b 2 because changing implies 
a change in bi or b 2 . Hence no clear direction in the shift of the predicted path can be observed by 
changing Hence no result is presented here. As will be seen later, knowledge of the effect of 
is not crucial to development of the calibration algorithms- the primary purpose of the thesis. 



Q X, Distance from left edge (m) 


Figure 3.1: Effect of a\ 
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Figure 3.4: Effect of a 4 


3.3 Effect of parameters of U 0 on predicted driver path 


The following are identified as parameters of U 0 : 04, m 0 , 64, 65, tp , t 1} and £2* The effect of these 
parameters on the predicted path are studied here. 

Effect of a 4 

Figure 3.4 shows the effect of a 4 on the predicted driver path on a 6 m road keeping all other 
parameters as constant. The value of the other parameters for the predicted paths in Figure 3.4 
are: (i) edge potential function parameters, a\ = 1.0, 02 = 4.0, bi = 0.32641, 62 = 0.04663 and (ii) 
obstacle potential function parameters, m* = 0.5, m a = 9.0, 64 = 2.0, 65 = 2.0, a = 0, /? = 0, 7 = 0, 
5 = 2.0, — 20, t 2 ~ 80 and tp = (2, 50). It can be seen that when <24 increases the deviation in the 

transverse direction (i.e. direction along the width) as well as deviation in the longitudinal direction 
(i.e. direction along the length ) increase. Further, the effect is symmetrical about a transverse axis 
thought the center of the obstacle. Also note that effect in the longitudinal direction is lesser than 
in the transverse direction. 
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Effect of m a 

Figure 3.5 shows the effect of m a on the predicted driver path on a 6 m road keeping all other 
parameters as constant. The values of the other parameters for the predicted paths in Figure 3.5 
are: (i) edge potential function parameters, a\ = 1.0, a 2 = 4.0, 61 = 0.32641, & 2 = 0.04663 and (ii) 
obstacle potential parameters, a 4 = 4.0 m* = 0.5, 64 = 2.0, 65 = 2.0, a = 0,/3 = 0, 7 = 0,5 = 2, 

= 20, t 2 = 80 and tp = (2,50). It can been seen that, when m a increases the transverse deviation 
(along the width of road) as well as deviation in the longitudinal (i.e. direction along the length) 
increase. However, the effect in longitudinal direction pronounced then in lateral direction; in fact, 
the maximum transverse deviation is insensitive to m a value. Further, the deviations are symmetric 
about a transverse axis through the center of the obstacle. 

Effect of 64 

Figure 3.6 shows the effect uf 64 on a 6 m road keeping all other parameters as constant. The 
values of the other parameters for the predicted paths in Figure 3.6 are: (i) edge potential function 
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parameters, a\ = 1.0, a 2 — 4.0, b\ = 0.032641, b 2 = 0.04663 and (ii) obstacle potential function 

parameters, a 4 = 2.0, m* = 0.5, m a = 8.0, ti = 10.0, 65 = 2.5, t 2 - 80.0 and tp = (50, 1.5). It can 

be seen that, when the value of 64 increases the transverse and longitudinal deviation on both side 
of tp decrease. The decrease, is however least at y = 50, the abscissa value for tp. 

Effect of 

Figure 3.7 shows the effect of & 5 on 6 m a road keeping all other parameters as constant. The 
values of the other parameters for the predicted paths in Figure 3.7 are: (i) edge potential function 
parameters, ai = 1.0, a 2 = 4.0, b\ — 0.032641, b 2 = 0.04663 and (ii) obstacle potential function 

parameters, a 4 = 2.0, mi = 0.5, m a = 8.0, b 4 = 2.50, ti = 20, t 2 ~ 75 and tp = (50, 1.5) It can be 

seen that, when 65 increases transverse and longitudinal deviation on both side of tp decrease. The 
decrease, is however least at y = 50, the abscissa value for tp. 

Effect of tp 

In order to study the effect of tp six cases are constructed. These are described in following table 
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x. Distance from left edge (m) 
Figure 3.7: Effect of 


Table 3.1: Cases for studying tp 


abscissa of tp 

ordinate of tp 

below 64 and 65 

above 64 and b§ 

in between 64 and 65 

above the y coordinate of obstacle 

Case I 

Case III 

Case V 

below the y coordinate of obstacle 

Case II 

Case IV 

Case VI 


The specific details of the cases and their corresponding figures are given below. 

Cases I and II 

Figure 3.8 shows five paths, each for a different value of the abscissa of tp. The obstacle is at (3, 50), 
hence, two of the paths are for values which place the tp above the y coordinate of obstacle (case 
I), two of the paths are for values which place the tp below the y coordinate of obstacle (case II), 
and one path is when the tp's abscissa value equal to the y coordinate of the obstacle. All the paths 
are plotted for the following parameter values: (i) edge potential function parameters, ai = 1.0, as 
= 4.0, b\ = 0.032641, 62 == 0.04663 and (ii) obstacle potential function parameters, a4 = 2.0, m l = 
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0.5, m a = 8 . 5 , 64 = 3.00, t x = 0.0, i > 5 = 3.0, t 2 = 100.0. It can be seen from figure that as tp moves 
away from the y coordinate of obstacle. The deviation reduce, irrespective of whether tp is above or 
below the position of the obstacle. 



x. Distance from left edge (m) 
Figure 3.8: Effect of tp in Cases I and II 


Cases III and IV 

Figure 3.9 shows the five paths, each for a different value of the abscissa of tp. The obstacle is at 
(3, 50), hence, two of the paths are for values which place the tp above the y coordinate of obstacle 
(Case III), two of the paths are for values which place the tp below the y coordinate of obstacle (case 
IV) and one path is w r hen the tp f s abscissa value equal to the y coordinate of the obstacle. All the 
paths are plotted for the following parameter values: (i) edge potential function parameters, a x = 
1.0, 0*2 = 4.0, b\ — 0.032641, b 2 = 0.04663 and (ii) obstacle potential function parameters, a 4 = 2.0, 
m t = 0.5, m a = 8.5, 64 = 3.00, t x = 0.0, 65 = 3.0, t 2 = 100.0. It can be seen from figure that as tp 
moves away from the y coordinate of obstacle. The deviation increase, irrespective of whether tp is 
above or below the position of the obstacle. 
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x. Distance from left edge (m) 
Figure 3.9: Effect of tp in Cases III and IV 


Cases V and VI 

In these Cases, the ordinate of tp is in between 64 and 65. For the same reasons as above five paths 
are plotted and showm in Figure 3.10. Each of these paths are plotted for the following parameter 
values: (i) edge potential function parameters, Oi = 1.0, a 2 = 4.0, bi = 0.32641, 62 = 0.04663 and 
(ii) obstacle potential parameters, G4 = 2.0, m* = 0.5, m a = 8.5, 64 = 3.50, t\ = 0.0, i> 5 = 1.5, 
£2 = 100.0. It can be observed from the Figure 3.10 that if the obstacle of tp is on the same side 
of the obstacle as the lesser of 64 and b$ (in this case 64 < 65) [notice here that 64 is the slope on 
the upstream side of obstacle and 65 is the slope on the downstream side of obstacle] then as tp 
moves away from obstacle the deviations increase. If, however tp is on the same side of obstacle as 
the greater of 64 and 65 (in this case 65) then as tp moves away from the obstacle the deviations 
decrease. 

Another aspect, which needs to be pointed out here is the effect of only changing the ordinate 
of tp keeping all other parameters fixed. Figure 3.11 shows four paths. Each is drawn for a different 
value of the ordinate. As can be seen, from the figure, the deviation decrease as the ordinate value 
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increases, irrespective of whether they are greater or smaller than 64 and 65. 

Effect of t\ 

Figure 3.12 shows the effect of t\ on the predicted path of driver on a 6 m road keeping all other 
parameters as constant. The values of the other parameters for the predicted paths in Figure 3.12 
are: (i) edge potential function parameters, a\ = 1.0, a<i = 4.0, £>1 = 0.32641, b 2 = 0.04663 and (ii) 
obstacle potential parameters, a 4 = 2.0, m* = 0.5, m Q = 8.5, 64 = 1.9, 65 = 1.9, t 2 = 85.0 and 
tp = (50,3). It can be seen that, when t\ moves away from the obstacle lateral and longitudinal 
deviations in the paths increase, although minimally. 



x, Distance from left edge (m) 
Figure 3.12: Effect of 


Effect of £2 

Figure 3.13 shows the effect of t 2 on the predicted path of driver on a 6 m road keeping all other 
parameter as constant. The values of the other parameters for the predicted paths in Figure 3.13 
are: (i) edge potential function parameters, a\ = 1.0, a 2 = 4.0, b\ = 0.32641, b 2 = 0.04663 and (ii) 
obstacle potential parameters, a 4 = 2.0, m* = 0.5, m a = 8.5, 64 = 1*90, t\ = 0.0, 65 = 1.9 and 
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tp = (50,3). It can be seen that, when £2 moves away from the obstacle deviations in the paths 
increase, although minimally. 



x. Distance from left edge (m) 
Figure 3.13: Effect of £2 


Based on the above parametric study, few important points relevant to the calibration process 
emerge. 

1. The parameters of the potential functions changes the predicted path in a variety of ways. 
Thus, it is felt that one can modify an initially predicted path to the observed path by varying 
these parameters. 

2. Further, it is possible to determine a direction of change of certain parameters given the 
observed and initially predicted path. For example, if the observed path lies to the left of the 
initially predicted path one can find that a\ has to be decrease, a<i has to be increase etc. 

3. There are however certain parameters for which it is not possible to predict any such direction 
of change. 

4. Lastly, one could identify certain parameters such as, £1 and £ 2 , which do not affect the predicted 
path much. 




Chapter 4 


PROPOSED CALIBRATION 
PROCESS 


In this chapter the proposed procedure for calibration are described. The chapter is divided into three 
sections in the first section overview of the proposed calibration procedure is explained. Section two 
describes Stage I of calibration procedure. Section three describes Stage II of calibration procedure. 


4.1 Overview of the proposed calibration procedure 


Given the observed path and the initial path (obtained using an initial set of values for the parameters 
of the potential field theory model) the calibration procedure uses the following steps to calibrate 
the parameter values of the potential field functions. 

1. Plot the initial path and observed path on a Cartesian space with abscissa, x, as the distance 
from the left edge of the road and ordinate, y, as the distance along the road from some 
arbitrary location. 

2. Divide the length of the road into a number of cross-sections (i.e. section with same ordinate). 
The cross sections can be arbitrarily close. 

3. For a cross section z, determine the error E{ as follows: 

E t = \(X' ob ,-X'J\ (4.1) 


IS 
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Figure 4.1: Stretch of the road illustrating the initially predicted path and the observed path 

where E x is the absolute error at i tk cross-section of the road, X l obs and X] n are the abscissa 
of observed and initial paths respectively at the i th cross-section. 

4. Check whether the E x at each cross-section is less than some threshold limit; if not, then 
correction in parameters of potential functions are required so go to next step; if yes, stop the 
process. 

5. Determine E t in the no-obstacle influence zone (see the Figure 4.1 for the definition of no- 
obstacle influence zone). 

6. Check whether the E % in the no-obstacle influence zone is less than the threshold limit; if 
no, apply Stage I of calibration procedure (as explained in Section 4.2); if yes, determine the 
characteristics of the observed path relative to the initial path after Stage I of calibration (as 
explained in Section 4.3) and apply Stage II correction rules to modify the parameter values 
of obstacle potential field function (as explained in Section 4.3.2). 

7. Replot initial path with modified (corrected) parameters of the potential functions and go to 
Step 2. These steps are shown schematically in Figure 4.2. 
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4.2 Calibration procedure for Stage I 

In Stage I, the corrections are applied to the parameters of the road edge potential functions to 
remove the errors in the no-obstacle influence zone (see the Figure 4.1 for the definition of no- 
obstacle influence zone) since in this zone only edge potential functions are active. The algorithm for 
Stage I corrections is based on the study performed in Chapter 3 on the effect of different parameters 
on the predicted path. As can be seen, from the study varying different parameters do have similar 
effects on the predicted path. For example, the predicted path shifts towards the right edge if a x is 
increased or a 2 is decreased, etc. Also, it can be observed that the rate of change of the predicted 
path with respect to different are different. Based on these observation, the algorithm is formulated 
thus: 

1. Determine X l obs and X\ n at i th cross-section of road in no-obstacle influence zone. 

2. Determine absolute error at i th cross-section in no-obstacle influence zone, I?*, as given below. 

E l = \{X i obs -Xl n )\ (4.2) 

3. Check whether E{ is less than threshold, if no, the correction in edge potential field function 
is required so go to next step; if yes, no correction is required so stop the process. 

4. If X l obs > X\ n increase a x and decrease a 2 (or if X l obs < X\ n decrease a x and increase a 2 ) 
with small step size without violating the following condition: a 2 > a i; and a 2 e~ b2W 4- a x < 
a x e~ blW -1- a 2 ). 

5. Either increase or decrease the value of b 2 with small step size (keeping in mind the conditions: 
a 2 > a x \ and a 2 e~ b2W 4- a x < a x e~ blW 4- a 2 ) so that error calculated by Equation 4.2 is further 
reduced and check whether error is less than threshold limit if yes, stop the procedure; if not, 
go to next step. Here, both increase and decrease of b 2 is tried because the position of the path 
also depends on and as such a clear understanding of the effect of b 2 alone is not available. 

6. Either increase or decrease the value of b x (keeping in mind the conditions: a 2 > a x ; and 
a 2 e~ b2W 4- a x < a x e~ blW 4- a 2 ) with small step size so that error calculated by Equation 4.2 
further reduce and check whether it is less than threshold limit if not, go to Step 4. Here, both 
increase and decrease of b x is tried because the position of the path also depends on and as 
such a clear understanding of the effect of b x alone is not available. 
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After applying Stage I corrections, in theory, the predicted path and the observed path in the no 
obstacle influence zone should become indistinguishable. Figure 4.3 shows a schematic of how the 
predicted would be in comparison to the observed path before and after Stage I calibrations. 




(a) (b) 

Figure 4.3: (a) Comparison of predicted path and observed path (a) before Stage I calibration and 
(b) after Stage I Calibration 


4.3 Calibration procedure for Stage II 

In stage two, corrections are applied to the parameters of obstacle potential field function. The 
calibration procedure can be divided into two parts: 

1. Measurement of characteristics of predicted path relative to observed path in obstacle influence 
zone after the completion of Stage I calibrations. 

2. Application of the rules of correction developed here to the parameters of the obstacle potential 


field function. 
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4.3.1 Measurement of characteristics of predicted path after Stage I rel- 
ative to observed path. 


The characteristics of the predicted path are measured in terms of the following factors: 

Initial path after edge correction 

Observed path 

5 in = Longitudinal deviation of initial path 

i 

x ob = Abscissa of observed path in section i 
x ' = Abscissa of intial path in section i 

in 

§ b = Longitudinal deviation of observed path 
t = Maximum transverse deviation of observed path 

max 



8. = X 
i obs 


m 


I e u = Maximum error in upper segement 

max 

£ = Maximum error in lower segment 

max 

P = Point of least error 


1 1.5 2 2.5 3 3.5 4 Total absolute error = £|( x ob S - x in )l 

x (along width) i 


Figure 4.4: Illustration of the characteristics of the observed path and predicted path after Stage I 
calibration 


1. Position of predicted path relative to observed path in a number of transverse sections in the 
obstacle influence zone. 

2. Total absolute error between observed and predicted path obtained using: 

Total absolute error = ]T Vt \{X' obs - X x in % where iL i” is a transverse section. 

3. Maximum transverse deviation of observed path and predicted path and their location along 
the road stretch consider(see Figure 4.4). 
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4. Longitudinal deviation of observed path and predicted path (see Figure 4.4). 

5. The point of least error (see Figure 4.4). 

6. Maximum error in the upper part of the path and maximum error in lower part of the path. 

These characteristics are measured in order help classify a particular scenario so that proper corrective 
actions can be taken. In the next section various scenarios (i.e. predicted path in comparison 
to observed path) are identified and the appropriate corrective actions in each of these cases are 
described. 

4.3.2 Correction rules for the parameters of potential field function of 
obstacles 

In this section the rules for the calibration of the parameter of the potential field function of static 
obstacle are presented. On studying the problem it was realized that a generalized calibration 
mechanism may not be feasible due to the complex nature of potential function and the way the 
parameters affect the path. Hence, it was decided to first classify a given condition (i.e., a given 
combination of the initially predicted path and the observed path) in to different “Situations.” The 
classification is based on the characteristics define in the earlier section. However, for clarity of 
understanding the different situation are presented here schematically in Figure 4.5. 

Classifications 

In all the parts of the figure the vehicle is assumed to travel from bottom to top of the figure. Part 
(a) of the figure is referred to as Situation 1. In this situation the initial path lies totally to the left 
of the observed path. That is, both the transverse and longitudinal deviation of initial path are more 
than that of the observed path. Part (b) of the figure is referred to as Situation 2. In this situation 
the initial path lies completely to the right of observed path. Part (c) of the figure is referred to 
as Situation 3. In this situation the initial path intersect the observed path. Further, the initial 
path is farther from the obstacle in the upstream side of the obstacle and closer to the obstacle in 
the downstream side of the obstacle when compared to the observed path. Part (d) of the figure 
is Situation 4 and represents the exact opposite of Situation 3. Part (e) of the figure represents 
Situation 5. In this*situation the initial path intersect the observed path twice. Further, the initial 
path is farther than the observed path in the vicinity of obstacle. Part (f) of the figure referred as 
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Figure 
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Situation 6 which is the opposite of Situation 5. 

It may be pointed out here that from the standpoint of determining the rules of the calibration, 
there exists no difference between the situations shown in Figure 4.5 and their mirror images about 
a longitudinal axis. In the next section the rules for calibration for each of these situation are 
developed. 

Rules of calibration 

Here, the rules of calibration for the different situations are described. These rules are obtained from 
the understanding of the effect of different parameters on predicted path (see Chapter 3). 

Situation 1: The rules are- (i) decrease a 4 , (ii) decrease m a , (iii) increase b 4) 65 and ordinate of 
tp such that the slope function f(y) takes a value greater than the initial value in every section, 
and (iv) if maximum error in upper segment, is more than that in lower segment, e l max , then 
decrease the abscissa value of tp and vice versa. Further, in this situation the rate of change of 
each parameter is kept same as them are equally effective in reducing in type of error seen in this 
situation. If however, e^ ax > e l max then 65 is change at a higher rate than b 4 and vice versa. 

Situation 2: The rules are- (i) increase a 4 , (ii) increase m a , (iii) decrease ft 4 , 65 and ordinate of tp 
such that the slope function f(y) takes a value smaller than the initial value in every section, and 
(iv) if maximum error in upper segment, is more than that in lower, e max , then decrease the 
abscissa value of tp and vice versa. Further, in this situation the rate of change of each parameter is 
kept same as them are equally effective in reducing in type of error seen in this situation, if however, 
e max > z l max then is change at a higher rate than b 4 and vice versa. 

Situation 3: Unlike in Situation 1 and 2 simple rules for changing the parameter like a 4 , m a , b 4: 
and tp are not possible in this case. The reason for this is that the initial path’s differences from the 
observed path have different characteristics in different segment. For example, the initial path has 
higher transverse deviation in the upstream side of obstacle but has lower transverse deviation in 
downstream side of obstacle when compare to observed path. Hence in this case the rules becomes 
slightly more complex. In following enumerated paragraphs the rules are described. 

(i) In every iteration change m a and a 4 in a manner such that the total absolute error (a 5 ex- 
plained in Section 5.14) reduces in that iteration. 
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(ii) Increase 64 and decrease 65 so that eventually the value of f(y ) is lower than initial value in 
the region where transverse deviation of initial path is lower than the observed path and vice versa. 

(iii) Change the abscissa value of tp such that total absolute error reduces. The ordinate of tp is 
changed based on two factors (i) which of |e^ aa: | and \c l max \ is greater, and (ii) whether the greater 
of the has a negative or a positive sign. If the greater of has a positive sign then reduce the ordinate 
tp , if however, it has a negative sign then the ordinate of tp is increased. 

Since parameter <24 and m a change the predicted path symmetrically with respect to transverse 
axis going through obstacle changes in these parameter play a small rule in reducing total absolute 
error. Hence in this situation parameter related to f(y) are changed at a higher rate than a 4 and 
m a . 

Situation 4 : Unlike in Situation 1 and 2 simple rules for changing parameter like a 4 , m a , 64, 65 
and tp are not possible in this case. The reason for this is that the initial path’s differences from the 
observed path have different characteristics in different segment. For example, the initial path has 
higher transverse deviation in the downstream side of obstacle but has lower transverse deviation 
in upstream side of obstacle when compare to observed path. Hence in this case the rules becomes 
slightly more complex. In following enumerated paragraphs the rules are described. 

(i) In every iteration change m a and a 4 in a manner such that total absolute error (as explained in 
section 5 . 14 ) reduces in that iteration. 

(ii) Decrease 64 and increase 65 so that eventually the value of f(y) is lower than the initial value 
in the region where transverse deviation of initial path is lower than the observed path and vice versa. 

(iii) Change the abscissa value of tp such that total absolute error reduces. The ordinate of tp 
is change based on two factors (i) which of |e“ aj .| and \e l max \ is greater, and (ii) whether the greater 
of the has a negative or a positive sign. If the greater has a positive sign then reduce the ordinate of 
tp, if however, it has a negative sign then ordinate of tp is increased. 

Since parameter a 4 and m a change the predicted path symmetrically with respect to transverse 
axis going through obstacle changes in these parameter play a small rule in reducing total absolute 
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error. Hence in this situation parameter related to /(y) are changed at a higher rate than a 4 and 
m a . 

Situation 5 : The rules are- (i) decrease a 4 (since maximum transverse deviation of initial path is 
more than of observed path), (ii) increase m a (because longitudinal deviation of initial path, 5 in is 
less than longitudinal deviation of observed path, 5 0 & s ), (iii) move the abscissa of tp in the direction 
of ordinate value of the point of maximum transverse deviation of the observed path, (iv) change i>4, 
65 and ordinate value of tp such that total error reduces. 

Since parameter m a is more effective in reducing the type of error seen in this situation. Hence 
in this situation m a is changed at higher rate than the other parameters. 

Situation 6: The rules are- (i) increase a 4 (since maximum transverse deviation of initial path is 
less than of observed path), (ii) decrease m a (because longitudinal deviation of initial, 6 in is more 
than longitudinal deviation of observed path, S 0 b s ), (iii) move the abscissa of tp in the direction of 
ordinate value of the point of maximum transverse deviation of observed path, (iv) change 64, 65 and 
ordinate value of tp such that total error reduces. 

Since parameter m a is more effective in reducing the type of error seen in this situation. Hence 
is this situation m a is changed at higher rate than the other parameters. 



Chapter 5 


RESULTS 


In this chapter results are presented to show that the proposed calibration mechanism can successfully 
calibrated the parameters of the potential functions. Thirteen different cases are studied here. 
Schematics of these cases are shown in Figure 5.1 as case (a) through case (m). The abscissa of 
each schematics is the distance along the width of the road from the left edge and the ordinate is 
the distance along the length of the road. The solid line in the figure is the predicted path from the 
potential field theory based model of driver behavior before calibration and the broken line is the 
observed path which can be obtained from real word observations. The solid object in the figure 
represents the obstacle. The motivation for choosing these cases for studying is that they represent 
all possible realistic combination of the initial predicted path and observed path. 

In the following, the results are shown for each of these thirteen cases. For each of these cases, 
expect for case (m) road width = 6m, with the left edge assumed at 0 m, obstacle is a circle of 
diameter 1 m and the center is placed at (50, 3). For case (m) only the placement of the obstacle is 
different; the center of the obstacle in this case is at (50, 2.3). The results for each case are presented 
using the following formate. First, a set of two figures marked (a) and (b) are presented. In the 
figures the solid line represents the path obtained using the potential field theoretical model with an 
initial (uncalibrated) set of parameter values. The dotted line is a hypothetical path representing (in 
theory) the observed path (which in practice will be obtained from real world data). The dashed line 
represents the path obtained using the potential field theoretical model with the calibrated parameter 
values. The figure marked (a) only presents the initial and observed paths. The figure marked (b) 
present the initial, observed and calibrated paths. 
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Second, for each case a table is provided which gives the initial parameter values and the calibrate 
values. 


1. Case (a) 

Figure 5.2 presents the result obtained using the proposed calibration mechanism. It c< 
be seen from figure that even though there exists reasonable difference between the initial ai 
observed path, after calibration the calibrated path is very close to the observed path. Table 5 
presents the initial and calibrated values of the parameters of the potential field theory mod< 


Table 5.1: Initial and calibrated parameter values for Case (a) 


Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 

a>\ 

1.0 

0.9110 

bi 

0.32641 

0.327332 

a 2 

4.0 

4.089002 

b 2 

0.04663 

0.045745 


2.0 

3.292 

m a 

8.0 

8.940 

b 4 

2.0 

1.952241 

tl 

0.0 

0.0 

a 

0.0 

-0.000005 

0 

0.0 

0.001139 

7 

0.0 

-0.071770 

<5 

2.0 

1.952241 

h 

2.0 

1.640931 

t'2 

100.0 

100.0 

tp (x, y) 

(50.0, 2.0) 

(42.10, 0.61) 


2. Case (b) 

Figure 5.3 presents the result obtained using the proposed calibration mechanism. It can 
seen from the figure that even though there exists reasonable difference between the initial a 
observed path, after calibration the calibrated path is very close to the observed path. Table i 
presents the initial and calibrated values of the parameters of the potential field theory mod 

3. Case (c) 

Figure 5.4 presents the result obtained using the proposed calibration mechanism. It can 
seen from the figure that even though there exists reasonable difference between the initial 





Figure 5.2: Result showing (a) the initial and observed paths and (b) the initial, observed and calibrated paths for Case (a) 
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Figure 5.4: Result showing (a) the initial and observed paths and (b) the initial, observed and calibrated paths for Case (c) 






Table 5.2: Initial and calibrated parameter values for Case (b) 
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Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 

ai 

1.0 

1.0890 

b i 

0.32641 

0.32553 

a 2 

4.0 

3.91100 

b 2 

0.04663 

0.047525 

<24 

2.0 

1.65 

m a 

8.0 

7.6500 

h 

2.0 

2.570301 

h 

0.0 

0.0 

a 

0.0 

-0.000004" 

p 

0.0 

0.000337 

7 

0.0 

0.002987 

s 

2.0 

2.570301 

h 

2.0 

2.274796 1 

t2 

100 

100.0 

tp (x, y ) 

(50.0, 2.0) 

(60.45, 3.10) 


and observed path, after calibration the calibrated path is very close to the observed path. 
Table 5,3 presents the initial and calibrated values of the parameters of the potential field 
theory model. 

Table 5.3: initial and calibrated parameter values for Case (c) 


Parameters of the 

Parameter values 

potential functions 

Initially assumed 

After calibration 

ai 

1.0 

0.9330 

bi 

0.32641 

0.32712 


4.0 

4.067002 j 

b 2 

0.04663 

0.045962 

0-4 

2.0 

1.3567 

ra a 

8.0 

6.3900 

b 4 

2.0 

3.729918 

1 1 

H 0.0 

0.0 1 

a 

0.0 

0.0 

p 

r 0.0 

0.000648 

7 

0.0 

-0.062219 

5 

2.0 

3.7299918 

&5 

2.0 

3.755372 


100.0 

100.0 

tp {x,y) 

(50.0, 2.0) 

(49.34, 2.21) 


4. Case (d) 

Figure 5.5 presents the result obtained using the proposed calibration mechanism. It can 
be seen from the figure that even though there exists reasonable difference between the initial 





Initial path 



Figure 5.5: Result showing (a) the initial and observed paths and (b) the initial, observed and calibrated paths for Case (d) 
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and observed path, after calibration the calibrated path is very close to the observed path. 
Table 5.4 presents the initial and calibrated values of the parameters of the potential field 
theory model. 


Table 5.4: Initial and calibrated parameter values for Case (d) 


Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 

ai 

1.0 

1.120 

b i 

0.32641 

0.325243 


4.0 

3.8800 


0.04663 

0.047835 

#4 

2.0 

3.6440 

m a 

8.0 

9.39992 

b 4 

2.0 

1.610029 

*i 

0.0 

0.0 

a 

0.0 

-0.000004 

e 

0.0 

0.000872 

7 

0.0 

-0.050558 

5 

2.0 

1.610029 

h 

2.0 

1.228277 

h 

100.0 

100.0 I 

*8“ 

li 

IT 

is 

(50.0, 2.0) 

(40.35, 0.723) 


5. Case (e) 

Figure 5.6 presents the result obtained using the proposed calibration mechanism. It can be 
seen from the figure that even though there exists reasonable difference between the initial and 
observed path, after calibration the calibrated path is very close to the observed path. Table 5.5 
presents the initial and calibrated values of the parameters of the potential field theory model. 

6. Case (f) 

Figure 5.7 presents the result obtained using the proposed calibration mechanism. It can be 
seen from the figure that even through there exists reasonable difference between the initial 
and observed path, after calibration the calibrated path is very close to the observed path. 
Table 5.6 presents the initial and calibrated values of the parameters of the potential field 
theory model. 

7. Case (g) 

Figure 5.8 presents the result obtained using the proposed calibration mechanism. It can be 
seen from the figure that even though there exists reasonable difference between the initial and 
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Table 5.5: Initial and calibrated parameter values for Case (e) 


Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 

ai 

1.0 

0.93300 

b i 

0.32641 

0.32712 

a 2 

4.0 

4.067002 

b- 2 

0.04663 

0.045962 

U4 

2.6 

1.867 

TTi a 

8.0 

8.950 

b 4 

2.0 

3.20001 

h 

0.0 

0.0 

a 

0.0 

0.000021 

P 

0.0 

-0.002864 

7 

0.0 

0.071480 

6 

2.0 

3.200001 

b 5 

2.0 

3.11968 

tp 

100.0 

100.0 

tp (x,y) 

(50.0, 2.0) 

(15.0, 3.7) 


Table 5.6: Initial and calibrated parameter values for Case (f) 


Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 


1.0 

1.120 

bx 

0.32641 

0.325243 

<*2 

4.0 

3.880 


0.04663 

0.047835 

(Z4 

2.0 

2.385 

m a 

8.0 

8.799910 

h 

2.0 

0.581535 

h 

0.6 

0.0 

a 

0.0 

0.000005 

7 

0.0 

0.000182 

p 

0.0 

-0.005492 

6 

2.0 

0.581535 

b 5 

2.0 

7.107835 

t2 

100.6 

100.0 

tp=(x,y) 

(50.0, 2.0) 

(10.40, 0.55) 





39 



(ui) qjgusj ppoy ‘A 




VO 

cn 


co 


8 

o 

W) 

x) 

<D 

4h 

<u 


vn 

h ge 


<N 


<D 

O 

0 


m 





Figure 5.6: Result showing (a) the initial and observed paths and (b) the initial, observed and calibrated paths for Case (e) 
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Figure 5.8: Result showing (a) the initial and observed paths and (b) the initial, observed and calibrated paths for Case (g) 
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observed path, after calibration the calibrated path is very close to the observed path. Table 5.7 
presents the initial and calibrated values of the parameters of the potential field theory model. 


Table 5.7: Initial and calibrated parameter values for Case (g) 


Parameters of the 

Parameter values 

potential functions 

Initially assumed 

After calibration 

CLi 

i.o 

1.120 

b i 

0.32641 

0.325243 


4.0 

3.8800 

£>2 

0.04663 

0.047835 

U4 

2.0 

1.88500 

m a 

8.0 

8.51 

b 4 

2.0 

9.162707 

*1 1 

0.0 

0.0 

a 

0.0 

-0.000012 

0 

o.o 

0.003222 

7 

0.0 

-0.289156 

<5 

2.0 

9.162707 

h 

2.0 

0.703910 

h 

100.0 

100.0 

75 

1! 

-& 

(50.0, 2.0) 

(85.10, 0.51) 


8. Case (h) 

Figure 5.9 presents the result obtained using the proposed calibration mechanism. It can be 
seen from the figure that even though there exists reasonable difference between the 

initial and observed path, after calibration the calibrated path is very close to the observed 
path. Table 5.8 presents the initial and calibrated values of the parameters of the potential 
field theory model. 

9. Case (i) 

Figure 5.10 presents the result obtained using the proposed calibration mechanism. It can be 
seen from the figure that even though there exists reasonable difference between the initial and 
observed path, after calibration the calibrated path 

is very close to the observed path. Table 5.9 presents the initial and calibrated values of the 
parameters of the potential field theory model. 

10. Case (j) 

Figure 5.11 presents the result obtained using the proposed calibration mechanism. It can be 
seen from the figure that even though there exists reasonable difference between the initial 
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Figure 5.10: Result showing (a) the initial and observed paths and (b) the initial, observed and calibrated paths for Case (i) 




Initial path J . _ L j Initial path 



Figure 5.11: Result showing (a) the initial and observed paths and (b) the initial, observed and calibrated paths for Case (j) 










Table 5.8: Initial and calibrated parameter values for Case (h) 


Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 

fli 

1.0 

0.9330 

b i 

0.32641 

0.32712 

a>2 

4.0 

4.067002 

h 

0.04663 

0.045965 

0,4 

2.0 

1.899 

m a 

8.0 

8.9500 

b 4 

2.0 

2.900171 

1 1 

0.0 

0.0 

a 

0.0 

-0.000027 

p 

0.0 

0.004457 

7 

0.0 

-0.180238 

<5 

2.0 

2.900171 

b 5 

2.0 

2.899993 

h 

100.0 

2.818590 

tp (x,y) 

(50.0, 2.0) 

(85.15, 3.40) 


Table 5.9: Initial and calibrated parameter values for Case (i) 


Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 

ai 

1.0 

0.95500 

b 1 

0.32641 

0.326892 


4.0 1 

4.045002 

62 

0.04663 

0.046185 

a.4 

2.0 

1.867 


8.0 

12.54 

b 4 

2.0 

1.30 

tl 

0.0 

0.0 

a 

0.0 

0.0 

p 

0.0 

-0.000517 

7 

" 0.0 

0.049223 

5 

^ 2.0 

1.300001 

b 5 

2.0 

1.30279 

t 2 

100.0 

100.0 

tp ( x,y ) 

(50.0, 2.0) 

(49.36, 2.5) 
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and observed path, after calibration the calibrated path is very close to the observed path. 
Table 5.10 presents the initial and calibrated values of the parameters of the potential field 
theory model. 


Table 5.10: Initial and calibrated parameter values for Case (j) 


Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 

a i 

1.0 

1.12 

h 

0.32641 

0.325243 

a 2 

4.0 

3.8800 

i>2 

0.04663 

0.047835 

O4 

2.0 

4.94300 

m a 

8.0 

5.09910 

h 

2.0 

2.79625 

h 

0.0 

0.0 

a 

0.0 

-0.000001 

v 

0.0 

0.001339 

7 

0.0 

-0.101512 

5 

2.0 

2.79625 

h 

1 2.0 

4.74830 

h 

1 100.0 

100.0 

tp (x,y) 

(50.0, 2.0) 

(40.21, 0.796) 


11. Case (k) 

Figure 5.12 presents the result obtained using the proposed calibration mechanism. It can be 
seen from the figure that even though there exists reasonable difference between the initial 
and observed path, after calibration the calibrated path is very close to the observed path. 
Table 5.11 presents the initial and calibrated values of the parameters of the potential field 
theory model. 

12. Case (1) 

Figure 5.13 presents the result obtained using the proposed calibration mechanism. It can be 
seen from the figure that even though there exists reasonable difference between the initial 
and observed path, after calibration the calibrated path is very close to the observed path. 
Table 5.12 presents the initial and calibrated values of the parameters of the potential field 

theory model. 
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Figure 5.13: Result showing (a) the initial and observed paths and (b) the initial, observed and calibrated paths for Case (1) 




Table 5. IT. Initial and calibrated parameter values for Case (k) 


Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 

a i 

1.0 

1.120 

b i 

0.32641 

0.325243 

a2 

4.0 

3.8800 

b 2 

0.04663 

0.427835 

<24 

2.0 

1.885 

m a 

8.0 

13.98600 

h 

2.0 

1.490010 

h 

0.0 

0.0 

a 

0.0 

0.000002 

0 

0.0 

-0.000610 

7 

0.0 

0.039099 

<5 

2.0 

1.490010 

b 5 

2.6 

1.361133 

*2 

100.0 

100.0 

tp = (x, 2/) 

(50.0, 2.6) 

(40.25, 2.21) 


Table 5.12: Initial and calibrated parameter values for Case (1) 


Parameters of the 
potential functions 

Parameter values 

Initially assumed 

After calibration 

Q>\ 

1.0 

0.9330 

b 1 

0.32641 | 

0.32712 

a 2 

4.0 

4.067002" 

62 

0.04663 

0.045962 

<24 

2.0 

4.0200 

m a 

8.0 

4.3900 

b 4 

2.0 

1.96097 

tl 

0.0 

0.0 

a 

0.0 

0.000001 

0 

0.0 

0.000372 

7 

*■“ 0.0 

-0.045790 

6 

2.0 

1.960979 

b 5 

2.0 

2.07192 1 

^2 

100.0 

100.0 

tp (x,y) 

(50.0, 2.0) 

(51.34, 0.721) 
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13. Case (m) 

Figure 5.14 presents the result obtained using the proposed calibration mechanism. It can be 
seen from the figure that even though there exists reasonable difference between the initial 
and observed path, after calibration the calibrated path is very close to the observed path. 
Table 5.13 presents the initial and calibrated values of the parameters of the potential field 
theory model. 


Table 5.13: Initial and calibrated parameter values for Case (m) 


Parameters of the 
potential functions 

Parameters values 

Initially assumed 

After calibration 

ai 

1.0 

2.8499 

h 

0.32641 

0.806805 


4.0 

3.45001 

b 2 

0.04663 

0.121544 

0-4 

2.0 

3.420 

m a 

8.0 

9.2100 

h 

2.0 

1.57306 

h 

0.0 

0.0 

a 

0.0 

0.0 

p 

0.0 

0.000216 

7 

0.0 

-0.024475 

<5 

2.0 

1.57306 

65 

2.0" 

• 1.56986 

h 

100.0 

100.0 

tp \x, y ) 

(50.0, 2.3) 

(51.43, 0.923) 


433-609 





(Ul) q}§U9{ pBO^I ‘A 



(Ul) 1{}§U9[ ptfO^ ‘/C 




Chapter 6 


CONCLUSIONS 


In this thesis a calibration mechanism which can calibrate the parameter values of potential field 
function of road edges and static obstacles is developed. In order to developed this mechanism a 
detail parametric study of potential field function is under taken and completed. Results showing 
the efficacy of proposed mechanism in calibrating the parameters are also presented. 

Although the proposed mechanism can ably calibrate the parameters it is felt that the efficiency 
of the algorithm can be further improved. For example, currently the step side for the corrections 
have been kept small. Although a attempt has been made to determine which parameter in which 
situation should be changed in larger step sizes, the determination is largely ad hoc. It is felt that 
a detail study on the step size should be conducted so that the calibrated values can be obtained 
faster. 


Another short coming of this thesis is that the observed paths were generated rather than 
experimentally obtained. Although it is felt that the proposed calibration algorithm will function 
equally well even if observed paths are experimentally obtained one can never be cent percent con- 
fident. 


It also felt the scope of this study should be enlarged to cover dynamic obstacles like moving 
vehicles, further, no attempt was made to calibrate the parameters of the models developed by 
Vasishta [2] and Agrawal [1] only the parameters of the potential functions were calibrated. This 
shortcoming must also be looked into. Also, the calibration was done based only on the observed 
position profile (i.e., the observed path of the driver) and did not use information on the velocity 
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profile (i.e., observed velocities of driver). This drawback needs to be rectified. 

•Finally a study should be conducted to see whether the calibrated potential field function can 
predict the driver path reasonably well for driving scenarios which are not used in the calibration 
process. This test is imperative to see whether the calibration process makes meaningful changes to 
the potential field functions. 
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